Implications of Simultaneous Requirements for Low Noise Exchange Gates in Double 

Quantum Dots 
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Achieving low-error, exchange-interaction operations in quantum dots for quantum computing 
imposes simultaneous requirements on the exchange energy's dependence on applied voltages. A 
double quantum dot (DQD) qubit, approximated with a quadratic potential, is solved using a full 
configuration interaction method. This method is more accurate than Heitler-London and Hund- 
Mulliken approaches and captures new and significant qualitative behavior. We show that multiple 
regimes can be found in which the exchange energy's dependence on the bias voltage between the 
dots is compatible with current quantum error correction codes and state-of-the-art electronics. 
Identifying such regimes may prove valuable for the construction and operation of quantum gates 
that are robust to charge fluctuations, particularly in the case of dynamically corrected gates. 



I. INTRODUCTION 

A quantum bit (qubit) typically encodes information 
in a two-level system. The exchange energy between 
quantum dots was first suggested as sufficient to per- 
form a universal gate set (unitary coherent manipulations 
of one and two qubits for logical operations) by Levy^i 
and subsequent exchange-based proposals for solid-state 
architectures have been suggested by Loss-DiVincenzo;^ 
Kane^ and Taylor 4 The exchange interaction causes a 
splitting between quantum states called the exchange en- 
ergy, which we denote J. Qubit rotations are performed 
experimentally by electrically increasing the exchange en- 
ergy for short times <^ 

Several important noise sources that can produce er- 
ror in the exchange operation include charge fluctua- 
tions {e.g. random telegraph noise, Johnson and shot 
noise), inaccuracy in electronics control (e.g. ringing and 
over/undershoot) and rotations due to inhomogeneous 
fields. 

Quantum error correction (QEC) schemes have been 
developed to cope with noise and errors in future quan- 
tum circuitsi^ A quantum error correction code intro- 
duces redundancy in the qubit information providing the 
ability to correct for errors through majority vote checks 
on the redundant basis bitsii These coding schemes are 
believed to be a necessary component of any future quan- 
tum computer because of the fragile nature of qubits. 
However, the codes provide benefit only for cases when 
the qubit gate error rate is less than a threshold value Pth 
above which the error correction circuit is more faulty 
than a bare qubit. Thresholds have been estimated for 
a number of cases and almost ubiquitously predict very 
strict limits on the tolerable error in the gate operations 
{e.g. Pth = 2x 10-^ and Pth = 2 x 10'^ from Refs.lH). 

A number of approaches are being pursued to minimize 
errors in qubit gates to achieve operations that are suffi- 
cient to realize the benefits of quantum error correction 
codes. The exchange gate couples the charge degree of 
freedom to the spin degree of freedom, which is useful for 
electronic control of the spins but also exposes the gate 



to errors induced by the electrodynamics of the system. 
A number of strategies have been proposed in the litera- 
ture to address different forms of errors {e.g. large J for 
fast rotations relative to noise sourcesji dJ/d{hia.s) « 
to suppress the impact of voltage fluctuations similar to 
those due to detuning^, and multiple rotation velocities 
for dynamically corrected gatingiS) . These strategies in- 
troduce a number of constraints, and it is not obvious 
that all of them can be simultaneously implemented in 
a double quantum dot given the physics of the system. 
In this paper we show that the simultaneous constraints 
arc consistent with a semi-qualitative model of a double 
quantum dot system. 



A configuration interaction (CI) method is used to 
study J as a function of parameters which specify a dou- 
ble quantum dot system. This CI method is more gen- 
eral than Heitler London (HL), Hund MuUiken (HM), 
and Hubbard model approaches, and is found invaluable 
to accurately calculate, in the single- valley case, energies 
for the bias range approaching and within the regime 
where there is two-electron occupation of one dot. The 
two-electron occupancy regime is relatively insensitive to 
the inter-dot bias {dJ /d{hias) « 0), making it an impor- 
tant regime to accurately calculate. Furthermore, this 
method is less computationally demanding than tech- 
niques requiring a large mesh, allowing a tractable search 
for robust exchange interaction parameters in the double 
dot system. 



We begin by describing our DQD model in section HIl 
and outlining the CI method used to solve it in section 
mil We then develop in section ITVl the constraints placed 
on a DQD's exchange energy by quantum error correction 
codes and controlling electronics. Results are presented 
in section |Vl and analyzed using the noise constraints. 
Finally, we discuss implications and a complementary ap- 
proach to noise mitigation in section IVIl and end with 
summary and conclusions in section fVIII 
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II. MODEL 

A lateral quantum dot singlet-triplet qubit qualita- 
tively similar to that described by Taylor et al<^ is ex- 
amined in this paper. To provide a semi-quantitative 
analysis we use gallium arsenide material constants m* = 
0.067 me and k = 12.9. The computational basis (i.e. the 
levels of the effective 2-state system) consists of the two- 
electron singlet and Sz ~ triplet states of lowest energy. 
J is the splitting between these two states. Note that 
the triplet states with = ±1 are split off with a dc 
magnetic field (typically of order 100 mT). The qubit 's 
effective many-body Hamiltonian is given by 

p2 e - ^ 

■H = Hi+H2 + ; m = - — + V{f) + — S-B (1) 

Kri2 2m* m* 

where ri2 = {fi — 7^2! and k is the GaAs dielectric con- 
stant. In "Hi, P = p — eA, andp = {px,Py) and f— {x,y) 
are the usual momentum and position operators of the 
i"^ electron. ^ is a vector potential for the magnetic field 
B = W X A, and V is the electrostatic potential. A con- 
stant perpendicular field B = Bz is considered here, and 
we restrict ourselves to two dimensions. 

The electrostatic potential V is generated by litho- 
graphically formed gates near the semiconductor inter- 
face. By applying different voltages to these gates at 
different times, the shape of V and the exchange energy 
can be tuned to perform operations on the qubit (e.g. see 
Ref.i. We idealize V as the minimum of two parabolic 
dots, 

V{x, y) = '^m*iJ- [min ((x - Bf + e, [x + Vf) + . 

^ .(2) 
The parameters e, L, and w, correspond to the bias, 
inter-dot distance, and frequency of the ground state as 
well as a measure of the confining well potential, respec- 
tively. Cuts of the 2D potential along the x-axis for dif- 
ferent e are shown in Fig. [T] One reason for choosing 
this parametrization is so that all three of these aspects 
can be varied independently. We also define the dot size 
d = ^Yi^ l{m*E[^ as the width of the ground state prob- 
ability distribution in a parabolic well with confinement 
energy Eq, where = hw. 

III. CALCULATION METHOD 

The CI method used to solve the system Hamiltonian 
^ can be decomposed into several steps. The initial and 
most complicated step is constructing a basis of uq s-type 
Gaussian functions. Each basis element g is parametrized 
by a position {xg,yg) and exponential decay coefficient Ug 

g{x, y) = iVe^"f (^-^f)'e-"9(?^~?^s)'e3#(y«^-^«f) (3) 

where A'^ is a normalization coefficient and B is the mag- 
netic field. A 2D mesh of points (given as input) speci- 
fies the positions of the elements in terms of two length 
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FIG. 1: DQD potential along x-axis, V{x,0), for L — 30 nm 
and Eq = 3.0 meV. The inset shows the potential obtained by 
solving Poisson's equation for both electrons and holes using a 
commercially available solver for an accurate DQD placement 
of gates and insulators. The dotted line is a double-parabolic 
potential and shows that a potential of this form is a relatively 
good approximation to the potentials expected in real devices. 



scales, ax and Uy, which can be thought of as effective 
lattice constants of the mesh. There is a mesh point at 
each of the dot centers, and we denote the exponential 
coefficient of the first basis element located at each of 
these points as ac- The exponential coefficient for the 
first element at all other points is denoted ao (generally 
«o 7^ etc)- When there are multiple basis elements at a 
point, the exponential coefficients of additional elements 
are found by multiplying the previous element's coeffi- 
cient by a constant factor A. Together the parameters 
Ux, dy, ac, ao, and A specify a Gaussian basis which is 
used in subsequent steps. The final values of these pa- 
rameters arc found by optimization of the full many-body 
energies, as described below. In all the results presented 
here — 36, with 18 elements on each dot arranged in 
two concentric 3x3 grids (so there are two elements on 
each dot center) , as shown in Fig. O 

Once a Gaussian basis is chosen, the single-particle 
Hamiltonian, minus the anomalous Zeeman term, is 
solved in this basis. The lowest n = 2nG (including 
spin degeneracy) of the resulting single-particle states 
are taken as an orthonormal (single-particle) basis, and 
all possible 2-particle Slater determinant states are con- 
structed, forming a riMB-dimensional two- particle basis, 
where umb = (2) ■ Lastly, the full many-body Hamilto- 
nian is diagonalized in this basis. This method consti- 
tutes a full CI with respect to the Gaussian basis. To 
improve convergence, the energies for the lowest singlet 
and Sz = triplet states are used to iteratively improve 
the Gaussian basis chosen initially. The final energies 
result from a Gaussian basis which is the direct sum of 
two smaller bases, one which minimizes the energy of the 
singlet and the other which minimizes the energy of the 
triplet. Complete details of this procedure will be given 
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FIG. 2: The relative positions of the Gaussian basis elements 
used in our CI method when no ~ 36. The center of each 
Gaussian is marked by a solid circle, and equipotential lines 
show the dot locations. The 18 elements on each dot are 
arranged in two concentric 3x3 grids which can have different 
spacings in the x- and j/-direction, as shown (there are two 
points at the center of each dot). The grid spacings are chosen 
to minimize the lowest singlet oi Sz = triplet energy, as 
described in the text. 



elsewhere 

Our results are only semi-quantitative, since the exact 
form of the potential is unknown and the problem is only 
approximately two-dimensional. However, the results 
give a more accurate qualitative and semi-quantitative 
picture than previous variational approaches, and are 
sufficient to resolve whether or not regimes exist which 
are robust to charge noise. Figure |3] below compares re- 
sults of our CI method with Heitler-London and Hund- 
Mulliken techniques. 



should therefore be targeted to make the ratio =j- as 

small as possible such that (■^)^ is at least smaller than 
the largest QEC threshold, that is, 

\AJ/Jf < Pth . (4) 

The magnitude of J is set by the gate time and target 
angle of rotation on the Bloch sphere. 



Use of the shortest gate times possible is a common 
strategy to minimize errors due to time-dependent de- 
coherence mechanisms. Electronics gate speeds are lim- 
ited by practical considerations {e.g. jitter and control of 
rise/fall times) and recent analysis suggests gate times in 
the range of 1 — 50ns<^i2ii^ The exchange energy must 
therefore be of order 1 — 0.02 ^eV for a 7r/2 rotation. 

B. No-op (idle) gate 

During a no-op gate, when no rotation is desired, J 
is ideally zero. If instead J takes finite value Jq, an er- 
roneous rotation ~ JaTa/h will occur over the time 
period tq. Inserting this into the condition |(5oP < Pth 
we find that the the magnitude of the exchange must 
satisfy 

\Jo? <n^Pthlrl: (6) 

where tq is length of time the qubit is idle. When Pth = 
10~^ this requires | Jo| < 6.5 — 0.13 neV for 1 — 50 ns gate 
times. 



C. Simultaneous Constraints 



IV. ERROR ANALYSIS 
A. Exchange (rotation) gate 

One model of an ideal exchange gate operation is to in- 
crease J from (near) zero to a finite value Ji for a time r, 
and then set J back to zero. This rotates the qubit about 
an angle 6 = Jir/ft, assuming that the exchange energy 
is dominant and therefore defines the axis of rotation. 

Constraints arise from the error thresholds demanded 
by error correction codes. In reality J cannot be perfectly 
controlled, and for an exchange error AJ the rotation 
angle becomes -\- 5, where 5 = AJt jfi. If J is intended 
to perform a 0-rotation, the gate time r = hO/ and 5 
is of order A J/ J. 

The probability of an error during an exchange gate 
can be estimated as approximately cos 5 5^ . The 
error probability should be engineered to be less than 
the predicted error threshold Pth of the quantum error 
correction, which was noted earlier to be dependent on 
the details of the QEC code and have a wide range of 
projected values {e.g. 10~^ — 10~^)i^i^ The quantity AJ 



Since an exchange gate operation involves tuning J 
between two values (usually zero and a finite value Ji ) , a 
robust gate requires the existence of at least two operating 
points where Eqs. H] and [5] are satisfied, respectively, for 
every needed rotation angle 6. 

Additional requirements for dynamically corrected gat- 
ing (DCG) are also potentially necessary in order to 
cancel other noise sources such as inhomogeneous quasi- 
static fields. At least three J(e) operating points are 
desired for DCG, and it is useful, but not necessary, for 
J to be negative at one of them.— A critical question 
theoretically is whether the dependence of J on the quan- 
tum dot properties can simultaneously realize several or 
all of these needs and thereby come closer to fulfilling the 
strict gate error requirements suggested by present error 
correction strategies. 

V. EXCHANGE ENERGY RESULTS 

The dependence of J on the electrostatic potential and 
magnetic field is examined using the CI method to idcn- 
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tify whether the predicted DQD exchange energy can 
meet the anticipated requirements for an exchange gate. 
We briefly consider the typical behavior of the exchange 
energy as a function of system parameters and then we 
analyze specific noise-robust regimes in detail. 



dJ/dB is larger for (0,2)-states, at large enough magnetic 
field the dip disappears. 



B. Noise- robust regimes 



A. Typical behavior 

Figure [3] shows the behavior of J as a function of e for 
varying perpendicular magnetic field strengths. When 
B < 1.1 T, the curve has two relatively flat sections at 
small and large e, where both the singlet and triplet states 
are in the (1,1) or (0,2) charge sector, respectively. When 
the inter-dot barrier is large enough, the quantum states 
of the DQD will often have an integral number of elec- 
trons in each dot. When there are ul and electrons 
in the left and right dots, respectively, the state is said to 
be in the [nL^nji) charge sector. Between the flat regions 
{e.g. in the case B = 1.1 T) J increases rapidly. This is 
because as e increases the singlet transitions to the (0,2) 
sector before the triplet, resulting in the DQD potential 
penalizing the triplet. The nearly constant value of J 
at large e is essentially the exchange energy of a doubly- 
occupied single dot with confinement energy Ef). 




FIG. 3: The exchange energy J as a function of inter-dot 
bias e for Eo = 3meV, L = 30 nm, and values of the mag- 
netic field around the point where the (0,2) singlet and triplet 
cross. At large and small e the curve J is relatively insensi- 
tive to changes in e. The transition from positive to negative 
J in the (0,2) regime (large e) gives rise to a local minimum 
near e = 3.5 meV. This provides a third region where J is 
relatively robust to e-variations. The curves labeled HL and 
HM show the Heitler London and Hund MuUiken methods for 
B = 1.1 T, respectively. 



Increasing the magnetic field favors the higher angular 
momentum of the triplet state relative to the singlet. The 
magnetic field needed to invert the singlet and triplet lev- 
els is lower for (l,l)-states than for (0,2)-states, and this 
creates a negative- J dip in the J vs. e curve at interme- 
diate magnetic fleld^' as shown in Fig. [31 Since the slope 



Given this general dependence of the exchange energy 
on inter-dot bias, three regimes can be identified which 
are relatively robust to e-noise: (I) at low-e, where the 
electrons are relatively isolated and J sa 0; (II) at high- 
e, where the singlet and triplet are in the (0,2) charge 
sector; and (III) at a local minimum, present for certain 
finite magnetic fields, where the singlet and triplet are 
between the (1,1) and (0,2) charge sectors. Whether the 
J vs. e curve in these regions can be made fiat enough 
to meet the stated requirements (Eq. 2] or [5]) is a central 
question we address in this work. 

Let us write the e-dependence of J explicitly, 
J — J(e), and define the average change in J(e) 
when e changes by Ae around eo: AJ(eo,Ae) ~ 
i (I J(eo + Ae) - J(eo)| + | J(eo) - ./(eo - Ae)|). Values 
of e used as rotation gate operating points must satisfy 
two criteria. Using Eq. [SJ 



|J(e)l < J,na.{e,T) = h\e\/T, 



(7) 



where, as previously, 9 is the gate rotation angle and r 
is gate time. Secondly, in order to satisfy Eq. |4l e must 
also be chosen so that 



Aecti < Aetargct = max { Ae 



AJ(e,Ae) 



J(e) 



< Pi 



th 



(8) 

where Aecti is the bias uncertainty achievable by the con- 
trolling electronics. Given r. 9, Pth, and e, Aetargct de- 
fines a target value of Ae which must be met in order 
to ensure noise robustness. In the case of a no-op gate, 
Eq. [S] needs to be satisfied with Jq equal to the exchange 
energy at and around the operating point, that is, 

max{|J(e')|' :e'e [e + Ae,e + Ae]} <ft2p,^/r2. (9) 

The quantity Aetargct can be defined in this case as the 
maximum value of Ae for which Eq. [2]is satisfied. Thus, 
for fixed Pth, Aecti, 9, and r, one simultaneously seeks 
large Aetargct and small \J\. 



1. Regimes I and II 

In regimes (I) and (II), Aetargct can be made arbitrarily 
large by increasing L or Eq, respectively. In both regimes, 
AJ/J decreases as the dots become more isolated, and 
eventually Eq. |S] will be satisfied. In regime (I), |J(e)| 
also decreases as the dots become more isolated, so that 
at large enough L oi Eq, Eq. [7] (or[9|) will be satisfied. In 
regime (II) , J takes a (usually non-zero) value dependent 
on Eq and B that can be made to satisfy Eq. [T] Thus, in 
regimes (I) and (II) one can theoretically satisfy Eqs. [7] 
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and [8l or Eq. [9] in the case of a no-op gate, by forming 
two isolated 1-electron dots or a single 2-electron dot. 

Consider the system given hy Eq ~ 3 meV, i = 60 nm, 
and B = 1667 mT. These parameters have been tuned 
within the range of physically reasonable values to result 
in an exchange curve with two relatively flat regions at 
low- and high-e. This curve is shown in Fig. 2] along 
with additional curves generated by varying Eq, L, and 
B around the point Eq = i meV, L = 60 nm, and B — 
1667 mT. 

At low-e (regime (I)), |J| < 10^^/ieV over a 3meV 
window around e = 0. Using Eq. [9] we find that the to- 
tal idle time tq can be up to (65.8 fis)\/Pth (66 ns when 
Pth = 10~^), which is typically the time of many gate 
operations. On the high-e flat (regime II), J w 1 ficV 
and dJ/de sa 5 x 10~^. This allows rotation gate times of 
order 1 ns, and from Eq.[5]Aetargot = (2 eV) y/Pth (2 nicV 
for Pth = 10~^, which practically is limited by the width 
of the flat region). Although these results are only semi- 
quantitative, this example shows that regimes exist for 
the DQD system where no-op and rotation operations are 
compatible with current quantum error correction and 
make realistic demands on controlling electronics tech- 
nology. This statement assumes, however, that only e- 
noise is present (i.e. parameters L, B, and Eq do not 
fluctuate). 

In actuality, L, B, and Eq cannot be perfectly con- 
trolled, and the variation of the exchange energy due 
to their fluctuations must be considered. In regime (I), 
when both the singlet and triplet are in the (1,1) charge 
sector, J is sensitive to the tunneling between the dots, 
and in general affected by L, Eq, and B. Because the ex- 
change energy is suppressed with increasing tunnel bar- 
rier, however, J and its variation over a given L-, Eq-, 
or _B-interval can be made arbitrarily small by choosing 
sufficiently large L, Eq and/or B. This is favorable for 
realizing a robust no-op. In the high-e regime (II), both 
electrons are almost completely confined to a single dot, 
and J is strongly dependent on Eq and B. The inter-dot 
spacing, on the other hand, has a relatively small effect on 
J that diminishes as e increases. In general, similar quali- 
tative behavior is obtained by increasing Eq or increasing 
B (due to their common confining effect on the electrons). 
This gives some freedom in selecting a dot size (Eq), and 
involves inherent trade-offs. For instance, in large dots 
smaller magnetic fields can accomplish the same effects, 
but larger dots are also more susceptible to disorder ef- 
fects {e.g. phonon induced spin-orbit couplingi^). 

The dependence of J on L, B, and Eq in regimes (I) 
and (II) can be seen in Fig. 0^, c, and e, respectively. 
Additionally, Fig. |4)d shows that the derivatives in these 
regimes are sensitive to L and therefore, even though 
changes in L do not affect the value of J on the upper 
flat, it must be sufficiently controlled that dJ/de remains 
within an acceptable range. Overall, to utilize the e- 
noise robustness of regime (II) requires an ability to hold 
L, B and Eq fixed to the extent that Eq. 2] is satisfied. 
Variations in L are least problematic, since keeping the 



dots sufficiently isolated will ensure that dJ/de is small. 
The typically strong linear dependence of J on and Eq , 
however, could not be avoided in the parameter ranges 
we studied. 

In summary, if L, B, and Eq can be held fixed precisely 
enough {i.e. to satisfy Eq. 2]), then there exist regimes 
of the DQD system which realize a robust no-op and 
rotation operation by varying only the inter-dot bias e. 
These robust regimes meet the requirements for current 
quantum error correction architectures and the control of 
e falls within current the capabilities of state-of-the-art 
electronics. 
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FIG. 4: Trends in the flat regions (regimes (I) and (II) in 
text) of exchange energy vs. bias curves around the point 
Eo = 3meV, L = 60 nm, and B = 1667 mT. Frames (a), 
(c), and (e) show the dependence on L (with B = 1667 mT 
and Eo = 3meV), B (with L = 60, and Eo = 3meV), and 
Eo (with L — 60 nm and B — 1667 mT), respectively. Frames 
(b), (d) and (f) show the derivative dJ/de of the curve at their 
left for e = and 15 meV. 



We note, however, that a large (0,2) singlet-triplet 
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splitting is necessary for loading and measurement so 
that the singlet can be selected with high probability. 
Thus, if the high-e regime is used for quantum opera- 
tions, there must be a method of temporarily increas- 
ing the (0,2) exchange splitting during initialization and 
measurement. 
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FIG. 5: Trends in relative curvature around the local mini- 
mum in exchange energy as a function of DQD bias e (regime 
II in text). Frames (a) and (b) sho-w the dependence on B 
(-with L = 30, and £^0 = 3 meV) and L (with B = 1300 mT 
and Eo = 3meV), respectively. Frames (c) and (d) plot the 
ratio A J/ J as a function of B and L, respectively, for the plots 
in frames (a) and (c) . It is seen that decreasing B or decreas- 
ing L correlates with smaller AJ/J, and this better robust- 
ness to e variations. In frame (e), curves with L — 60 nm and 
similar values at their minima, but different Eo, are shown, 
and it is seen that decreasing Eo (and B) lead to greater e- 
robustness. 



2. Regime III 

In regime (III), dJ/de ~ exactly at the minimum, 
and it more informative to study the behavior of Aetargct 
at the minimum, "which is a measure of the curvature of 
J(e) relative to its magnitude. For a given Pth, Aetargct 
can be increased by decreasing B, L, or Eq. Decreas- 
ing B gives less energetic advantage to the triplet state 
(relative to the singlet) and results in the minimum be- 
coming less sharp as well as occurring at smaller |J|, as 
shown in Fig. [SK- Overall, Aetarget increases as seen in 
Fig. [5b • Either decreasing L or Eq increases the overlap 
between (1,1)- and (0,2)-states. Even though this pushes 
the minimum to larger | J| because the (1,1) states have 
greater (negative) exchange energy, Aetargct increases be- 
cause the (l,l)-(0,2) transitions of the singlet and triplet 
occur more gradually, and are farther separated in e. 
This dependence on L is shown in Fig. [SJ; and d, and 
the dependence on Eq in Figs. [5^ and f. Note that in 
Fig. [5] we compensate an increase in Eq by reducing B 
to obtain local minima which occur at similar values of 
J. Even without this compensation, lower values of Eq 
give larger Aetargct for fixed Pth- Because the singlet 
and triplet states have mixed (1,1) and (0,2) character in 
regime (III), variations in L, Eq, and B are particularly 
effective at changing the charge distribution of the sin- 
glet and triplet, and thereby the exchange energy. The 
sensitivity of J to variation in L is greater in this region 
than in either the low- or high-e region, and the sensitiv- 
ity to Eq and B lies between that of the low- and high-e 
regions. 

Ideally, all three regimes could be used as robust op- 
erating points simultaneously with the low-e region serv- 
ing as an idle point. Multiple rotation speeds and nega- 
tive rotation are potentially desirable for DCG. Achieving 
this goal is challenging for several reasons. The first is 
that broadening the regime- (III) minimum by coupling 
the dots more strongly is correlated with the exchange 
energy at e = increasing in magnitude (see Fig. [S]). 
This sets up a competition between a robust idle state 
in regime (I) requiring small \J\, and a robust operating 
point in regime (III), which requires a broad minimum. 
Indeed, we find that for Pth = 10^^ and r ~ 1 ns it 
is impossible to satisfy the constraints for an idle gate in 
regime (I) and a rotation gate in regime (III) at the same 
time. Secondly, achieving a given value of \J\ at large-e 
requires tuning either the dot size (via Eq) or magnetic 
field. But since the magnetic field must be tuned to give 
a usable minimum in regime (III), tuning the dot size 
is necessary to simultaneously achieve a usable regime 
(II). For |J| of order /icV these constraints can lead to 
very large dots {d greater than 200 nm). Figure [H] shows 
an example of a curve containing operation points in all 
three regimes (and with regime I an idle point). Utilizing 
subsequent error correction, however, requires very fast 
gating of 0.1 ns and high detuning accuracy of 0.02 meV, 
which represents a significant technological challenge. In 
the end, we find that while the physics of the DQD system 
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allows at least three robust operating points reachable by 
only changing the inter-dot bias, utilizing all three will 
require either finer electronics control or more efficient 
quantum error correction algorithms. 




FIG. 6: Exchange energy J vs. bias e with tliree potential op- 
erating points lying at e ~ 0, 1.6, and 5meV. To be compati- 
ble with an error threshold of 10"*, however, would require a 
gate time of order 0.1 ns and control of e to a few hundredths 
of a micro-electronvolt. 



VI. DISCUSSION 

It may not be necessary, however, to have three (or 
more) operating points separated only by changes in the 
inter-dot bias. One could envision implementing two ro- 
tation speeds in a DQD qubit by varying another pa- 
rameter, such as Eq or L, as well as e, and to perform 
this alternate variation while e = and the qubit is in a 
robust no-op state. Or perhaps only one rotation speed 
will be necessary to begin with. Our results indicate that 
a robust no-op should be accessible using current control 
and error correction technology, and that a robust ro- 
tation operation is also feasible as long as the shape of 
and spacing between the dots can be controlled with high 
precision. 

Though dynamics are not studied in this work, it is im- 
portant to realize that utilizing a (0,2)-flat region (regime 
II) requires the inter-dot bias to be quickly changed so 
that relatively little time is spent in the region of the 
J curve between the low- and high-e flats. The speed 
at which this bias change occurs is limited by the gap 
to higher energy singlet and unpolarized triplet states as 
dictated by the adiabatic theorem. We have considered 
such restrictions, and find that the gap to excited levels 
remains large enough that the qubit can be moved adia- 
batically between regimes I and II with the vast majority 
of the gate time spent in the noise robust regimes. This 
docs, however, set a bound on how weakly the dots can 
be coupled, since the gap to excited states decreases with 
the inter-dot coupling. 

The architecture of the DQD can be used to mitigate 
the effects of the exchange energy's sensitivity to L and 
£^0- To the extent that the actual DQD potential re- 
mains a double-parabolic well, there will be a mapping 



from sets of gate voltages to the parameters e, L, and 
Eq. Variations in these parameters is thus determined 
by their dependence on the gate voltages which vary to 
perform a qubit operation. In this work, we have identi- 
fied regions of (e,L,£'o)-space which are favorable for sup- 
pressing charge noise because they are flat, or nearly flat, 
along at least the e-direction. By modifying the architec- 
ture of a DQD device, one can hope to map the pathways 
in gate- voltage space that perform qubit operations onto 
pathways in (e,L,£^o)-space that begin and end along flat 
regions. In the present work we specifically focus on ro- 
bustness to e variations, and the ideal architecture would 
allow gate voltages to change e while keeping Eq and L 
fixed. For example, to keep Eq fixed for the right dot, 
changes in inter-dot bias might be controlled exclusively 
by varying the voltages of gates around the left dot. Such 
architecture engineering was first proposed by Friesen et 
al.ji^ where by moving the electrons along parallel chan- 
nels instead of directly toward or away from each other, 
the inter-dot separation L varies only quadratically in 
the gate voltages at a finite- J operating point, instead of 
linearly. In that work, however, since the DQD is in a 
regime where J depended exponentially on L, the archi- 
tecture serves only to reduce the sensitivity of L to the 
gate voltages, not to map the pathway onto a flat curve 
in (e^LjEa) parameter space. 



VII. CONCLUSION 

In summary, configuration interaction calculations on 
a singlet-triplet DQD qubit using GaAs material parame- 
ters have been carried out. Three regimes have been iden- 
tified in which the exchange energy is relatively insensi- 
tive to changes in the inter-dot bias e. The CI method is 
necessary to both qualitatively and quantitatively calcu- 
late the dependence of J on critical parameters such as 
detuning e, dot energy Eq, and dot separation L, com- 
pared to previous more approximate schemes such as HL 
or HM. In particular, the CI method is found invalu- 
able for calculations of critical regions such as when the 
dots are strongly coupled or when a single dot is doubly- 
occupied. Namely, it captures the regime in which both 
the singlet and triplet transition into the (0,2) charge 
sector. 

By tuning only the intcr-dot bias it is possible to travel 
between two or possibly three (with advances in electron- 
ics technology) of these robust regimes, which is desired 
for dynamically corrected gates and suggests how they 
might be implemented. These types of calculations are 
needed to provide guidance regarding accuracy require- 
ments for e, Eq, and L given a QEC threshold. We note 
the adverse effects caused by the sensitivity to certain 
parameters may be avoided by clever design of the qubit 
control electronics and architecture. 
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